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A STRATEGY OF NUMERIC SEARCH FOR 
PERFECT CUBOIDS IN THE CASE OF 
THE SECOND CUBOID CONJECTURE. 


A. A. Masharov, R. A. Sharipov 


Abstract. A perfect cuboid is a rectangular parallelepiped whose edges, whose face 
diagonals, and whose space diagonal are of integer lengths. The problem of finding 
such cuboids or proving their non-existence is not solved thus far. The second cuboid 
conjecture specifies a subclass of perfect cuboids described by one Diophantine equa¬ 
tion of tenth degree and claims their non-existence within this subclass. Regardless 
of proving or disproving this conjecture in the present paper the Diophantine equa¬ 
tion associated with it is studied and is used in order to build an optimized strategy 
of computer-assisted search for perfect cuboids within the subclass covered by the 
second cuboid conjecture. 


1. Introduction. 

For the history and various approaches to the problem of perfect cuboids the 
reader is referred to [U42]. In this paper we resume the research initiated in [43-47]. 
The papers [48-60] deal with another approach based on so-called multisymmetric 
polynomials. In this paper we do not touch this approach. 

Perfect cuboids are described by six Diophantine equations. These equations are 
immediate from the Pythagorean theorem: 


x\+xl+xl- =Q, xl + xl-dl=Q, 

x\ + x\ — d 2 = 0 , x\ + X 2 — d^ = 0 . 

The variables Xi, X 2 , X 3 in (l.I) stand for three edges of a cuboid, the variables di, 
d 2 , ds correspond to its face diagonals, and L represents its space diagonal. 

In [43] an algebraic parametrization for the Diophantine equations (1.1) was 
suggested. It uses four rational variables a, /3, v, and z: 


xi 2v 

T ^ TT^’ 


di 1 — v'^ 

T ^ 1+^2’ 


X2 2z{l — ^ 2 ) 

T ^ (l+u2)(l-f z2)’ 


X3 _ (1 - (1 - Z^) 

L ~ (l + ^; 2 )(l + z 2 )’ 


d 2 (1+u2)(1+z2)+2z(1-u2) 

L (l+^; 2 )(l+^ 2 ) P’ 


ds _ 2 (u2 2:2 _|_ 

L (1-I-u2) (1 _|_ 2:2) 
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The variables a, j3^ v in (1.2) are different from the original ones which are used in 
[43], here we use a and /3 instead of a and &, and we use v instead of u. 

Only two of the four variables a, (3, v, and 2 are independent. The variables a 
and j3 are taken for independent ones. Then the variable v is expressed through a 
and /? as a solution of the following algebraic equation: 


/3^ + (6 - 2 - 2 a" /3^) + (4 /3^ a^ + 

+ 4 — 12/3^ + /3^ + /3^) + (6 + 6 v'^ f3‘^— 

- 8 ;32 ^2 _ 2 ^4 ^2 _ 2 ^4 ^2 _ 2 q ,4 ^2 _ 2 ^4 ^ 2 ) 

+ a^ + 4a2^;2 + 4^2^2 _ ^2^32 ^2^ (5 ^2 _ 2 ^2 _ 2 ^2^ ^ ^ 


(1.3) 


Once the variable v is expressed as a function of a and j3 by solving the equation 
(1.3), the variable z is given by the formula 


(1 + v'^) (1 - /3^) (1 + Qf^) 


(1.4) 


2 (1 +/?2) (1 — q; 2 .y2^ 

The equation (1.3), along with the formula (1.4), produces two algebraic functions 
v = v{a,/3), z = z{a,l3). (1.5) 

Substituting (1.5) into (1.2), we get six algebraic functions 


xi = xi{a, 13, L), 
X2 = X2ia,/3,L), 
X3 = X3{a, 13, L), 


di = di{a,f3,L), 
d2 = d2(a,/3,L), 
ds = d3{a,j3,L), 


( 1 . 6 ) 


which are linear with respect to L. The functions (1.6) satisfy the cuboid equations 
(1.1) identically with respect to a, (3, and L. This fact is presented by the following 
theorem (see Theorem 5.2 in [43]). 

Theorem 1.1. A perfect cuboid does exist if and only if there are three rational 
numbers a, jd, and v satisfying the equation (1.3) and obeying four inequalities 
0<a<l, 0</3<l, 0<t><l, and {a + 1) (/3 + 1) > 2. 

The rational numbers a, (3, and v cam be brought to a common denominator: 


a = 


P = 


V = 


t t t 

Substituting (1.7) into (1.3), one easily derives the Diophantine equation 

+ (6 — 2 — 2 6^) + (u^ + 6^ + + 4 + 

+ 4 52 y2 _ 42 ^2 ^2) ^8 (6 y2 g ^2 ^4 _ g ^2 ^2 ^ 2 _ 

— 2 — 2 6^ — 2 — 2 5^ a^) + (4 5^ a^+ 

+ 4:a'^u^b^ — 12 b^ + u'^ a'^ + u'^ 6^ + b'^) 

+ (6 u ^ b ^-2 a ^ b ‘^-2 6 ^) = 0 . 


(1.7) 


( 1 . 8 ) 


Theorem 1.1 then is reformulated in the following form (see Theorem 4.1 in [44]). 



A STRATEGY OF NUMERIC SEARCH FOR PERFECT CUBOIDS ... 


3 


Theorem 1.2. A perfect cuboid does exist if and only if for some positive coprime 
integer numbers a, b, and u the Diophantine equation (1-8) has a positive solution 
t obeying the inequalities t > a, t > b, t > u, and {a + t) (b + t) > 2t^. 

In [44] the Diophantine equation was treaded as a polynomial equation for t, 
while a, b, and u were considered as parameters. As a result in [44] several special 
cases of the equation (1.8) were specified. They are introduced through the following 
relationships for the parameters a, 6 , and u: 

1) a = b^u] 3) bu = a^\ 5) a = u ^ b; 

(1.9) 

2) a = b = u; 4) au = b^\ 6 ) b = u ^ a. 

The cases 2, 5, and 6 are trivial. They produce no perfect cuboids (see [44]). The 
case 1 corresponds to the first cuboid conjecture (see [44]). It is less trivial, but it 
produces no perfect cuboids either (see [45]). The cases 2 and 4 correspond to the 
second cuboid conjecture (see [44] and [46]). The case, where none of the conditions 
(1.9) is fulfilled, corresponds to the third cuboid conjecture (see [44] and [47]). 

In this paper we consider the cases 3 and 4 associated with the second cuboid 
conjecture. In the case 3 the equality bu = of is resolved by substituting 

a=pq, b = p^, u = q^. (1-10) 

Here p ^ q are two positive coprime integers. Upon substituting (1.10) into the 
equation ( 1 . 8 ) it reduces to the equation 

{t - a) (t + a) Qpq{t) = 0 ( 1 - 11 ) 

(see [45]), where Qpq{t) is the following polynomial of tenth degree: 

Qpq(t) = t^° + (2g^ +P^) [iq^ - 2p^)t® + (g® + 10p^g®+ 

+ 4p^g'‘ - 14p®g2 _^p®)t® _ 14^2 ^6 4^4 ^4 42 ) 

+ 10 P® p 8 ) ^4 _ p 6 ^6 (-^2 2 p 2 ) ( 3 p 2 _ 2 q^) f _ glOplO, 

The case 4 is similar. In this case the equality au = b^ is resolved by substituting 

a = p^, b = pq, u = q^. (1-13) 

Upon substituting (1.13) into the equation (1.8) it reduces to the equation 

{t-b){t + b)Qpq{t) =0. (1-14) 

The roots t = a, t = —a, t = b, and t = —6 of the equations (1.11) and (1.14) do 
not produce perfect cuboids (see Theorem 1.2). Upon splitting off the linear factors 
from (1.11) and (1.14) we get the equation 

t^° + ( 2 g 2 +p 2 ^ ( 3^2 _ 2 p 2^^8 (^8 10 ^ 2^6 _^4p4 ^4 _ 

- 14p®g2 _^p®)t® -p 2 g 2 (^8 _ I 4 p 2^6 _^4p4^4 I 0 p®g 2 + 

+p®) t'^ - p® g® {q^ + 2p2) (3p2 - 2 q^) f - q^°p^° = 0. 


(1.15) 
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Conjecture 1.1. For any positive coprime integers p ^ q the polynomial Qpq{t) 
in (1-12) is irreducible in the ring 1[t]. 

Conjecture 1.1 is known as the second cuboid conjecture. It was formulated in 
[44]. In particular it claims that the equation (1.15) has no integer roots for any 
positive coprime integers p ^ q. We do not try to prove or disprove Conjecture 1.1 
in this paper. Instead, we study real positive roots of the equation (1.15) in the 
case where q is much larger than p. Using asymptotic expansions for the roots of 
the equation (1.15) as g —>■ +oo, below we build an optimized strategy of computer- 
assisted search for perfect cuboids in the realm of Conjecture 1.1. 

2 . Asymptotic expansions for roots of the polynomial equation. 

Note that the polynomial Qpq(t) in (1.12) is even. Along with each root t it has 
the opposite root —t. We use the condition 

f t > 0 if t is a real root, 

{ ( 2 - 1 ) 

( Re(t) ^ 0 and Im(t) >0 if t is a complex root, 

in order to divide the roots of the equation (1.15) into two groups. We denote 
through ti, t 2 , t 3 , ti, ts the roots that obey the conditions (2.1). Then te, h, t%, 
tg, tio are opposite roots of the equation (1.15): 


tg — —ti, t-j — —tg, ts — —ta, tg — —t 4 , tiQ — —t^. (2-2) 


Typically, asymptotic expansions for roots of a polynomial equation look like 
power series (see [61]). In our case we have the expansions 


U{p,q)=aq^* 


1 + 




as g —>• -\-oo. 


(2.3) 


The coefficient Ci in (2.3) should be nonzero: Ci ^ 0. 

Let’s substitute (2.3) into the equation (1.15). For this purpose we represent the 
polynomial Qpq{t) from (1.12) formally as the sum 


10 10 

Qp,{t)=Y.T.^^r{p)q^t^. (2.4) 

m—0 r—0 


Each nonzero term in (2.4), i. e. a term Amr(p) q^ with the nonzero coefficient 


yields the sum 


Amr(p) 0; 


5^,(p,g) = A^,(p) (2.5) 

s<m ai+r 

Taking into account (2.5), the equation (1.15) is written as 

10 10 

E '^^^r{p,q) = 0. 

m—0 r—0 


( 2 . 6 ) 
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The equality (2.6) should be fulfilled identically with respect to the variable 
q —>■ + 00 . Since Ci ^ 0, a necessary condition for that is the coincidence of 
exponents of at least two summands of the form Smr{q) in the leading order with 
respect to the variable q. This yields the equalities 

mi a* + ri = m 2 ai+r2 = s^ax- (2.7) 

The maximality of the exponent in (2.7) means that all exponents are not greater 
than Smax, i-e. we have the following inequality: 

moi + r^s^ax for all r and m such that Amrip) Q- (2-8) 

Lemma 2.1. The coincidence mi = m 2 in the formula (2.7) is impossible. 

Proof. Indeed, due to (2.8) the coincidence mi = m 2 would mean ri = r 2 . But the 
sum (2.6) has no two summands with simultaneously coinciding indices r and m. 
Lemma 2.1 is proved. □ 

Let’s treat m and r as coordinates of a point on the coordinate plane. Since m 
and r are integer, such a point belongs to the integer grid, being its node. The 



Fig. 2.1 


numbers mi, ri and m 2 , r 2 from (2.7) mark two nodes of this grid. These are the 
points A and B in Fig. 2.1. Due to Lemma 2.1 from the equality (2.7) we derive 
the following formula for the exponent ap. 


r2 - ri 

(Xi — 

Vfl2 — TUi 


(2.9) 


The right hand side of the formula (2.9) up to the sign coincides with the slope of 
the straight line connecting the nodes A and B in Fig. 2.1: 


Oii — kAB- 


(2.10) 
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The nodes A and B correspond to some nonzero summands in the sum (2.4) 
being a formal presentation of the polynomial (1-12). They are selected by the 
maximality condition for the parameter s = mai + r. The maximum is taken over 
all summands in the sum (2.4) for a fixed value of ai. 

Lemma 2.2. The exponent in the asymptotic expansion (2.3) is determined by 
the slope of a straight line connecting some two nodes of the integer grid associated 
with some two nonzero terms in the polynomial (2.4). 

Let C be some node of the integer grid in Fig. 2.1 associated with some nonzero 
summand of the sum (2.4) and different from the nodes A and B. Its coordinates 
m and r satisfy the inequality (2.8). From (2.7) and (2.8) one derives the inequality 


m ai + r ^ mi Oi + ri. 


Let’s write this inequality as follows: 


Oj (m - mi) ^-(r - ri). (2-11) 

In Fig. 3.1 three positions of the node C relative to the node A are shown. The 
node C can be located to the left of the node A, to the right of the node A, or on 
the same vertical line with the node A. In the first case m < mi. In the second 
case m > rrii. And finally, in the third case m = mi. 

In the first case, i.e. if m < mi, from (2.11) we derive 

r — 7*1 

a, ^-(2.12) 

m — mi 

The right hand side of the inequality (2.12) up to the sign coincides with the slope 
of the line AC. Applying (2.10), we get the inequality —kAB ^ —kAC- Inverting 
signs, we write this inequality in the following form: 

kAC ^ kAB- (2-13) 


In the second case, i. e. if m > mi, from (2.11) we derive 

a, < (2.14) 

m — mi 

By analogy with (2.13) the inequality (2.14) is transformed to 


kAC ^ kAB- (2-15) 

And finally, in the third case, i. e. if m = mi, from the inequality (2.11) we derive 

0<-(r-ri). (2.16) 

The inequality (2.16) is equivalent to the following inequality: 


r ^ ri. 


( 2 . 17 ) 
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Each of the inequalities (2.13), (2.15), and (2.17) in its case means that the point 
C is located not above the line AB. This fact is formulated as a lemma. 

Lemma 2.3. All nodes {m,r) of the integer grid associated with nonzero sum¬ 
mands in the polynomial (2.4) are located not above the line AB on which the nodes 
implementing the maximum of the parameter s = m + r are located. 

In order to apply Lemmas 2.1, 2.2, and 2.3 let’s mark all of the nodes associated 
with the polynomial (1.12) on the coordinate plane. 



Fig. 2.2 


Definition 2.1. For any polynomial of two variables P{t, q) the convex hull of all 
integer nodes on the coordinate plane associated with monomials of this polynomial 
is called the Newton polygon of P{t,q). 

Remark. Note that in our case the polynomial (1.12) depend on three variables 
p, q, and t. However, we treat p as a parameter and consider Qpqit) as a polynomial 
of two variables when applying Definition 2.1 to it. 

The Newton polygon of the polynomial (1.12) is shown in Fig. 2.2. Its boundary 
consists of two parts — the upper part and the lower part. The upper parts is 
drawn in green, the lower part is drawn in red. In Fig. 2.2 the nodes on the upper 
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boundary of the Newton polygon are denoted according to the formula (2.4). The 
coefficients Amr{p) in (1-12) associated with these nodes are given by the formulas 

4loiO = —yl210 = 2p®, ^410 = — 

(2.18) 

^68 = 1, ^84 = 6, Aioo = l- 

Theorem 2.1. The values of exponents in the expansion (2.3) for roots of the 
equation (1.15) are determined according to the formula at = —k, where k stands 
for slopes of segments of the polygonal line being the upper boundary of the Newton 
polygon in Fig. 2.2. 

Theorem 2.1 is immediate from Lemmas 2.2 and 2.3. The formula = — fc in 
this theorem follows from the formula (2.10). In our particular case we have 


a* =0, a, = 1, tti = 2. (2-19) 

The options (2.19) are derived from Fig. 2.2 due to the above theorem. 

3 . Leading terms in asymptotic expansions. 

The term Ci ( 7 “* obtained upon expanding brackets in (2.3) is called the leading 
term of the asymptotic expansion (2.3). Three options for the value of ai are given 
by the formula (2.19). Let’s consider each of these options separately. 

The case ai = 0. This case corresponds to the horizontal segment on the upper 
boundary of the Newton polygon in Fig. 2.2. This segment comprises three nodes 
4 I 410 , ^ 210 , and ^0 10 - Therefore, substituting the expansion (2.3) with 0 ^ = 0 
into the equation (1.15), we get the following equation for Cp. 

A 4 10 + A 2 10 + Aq 10 = 0. (3-1) 

Taking into account (2.18), the equation (3.1) is transformed to 

p'^a^ -2p^a^ +p^° = 0. (3.2) 

The equation (3.2) has two real roots 




C. = -p^ 


(3.3) 


each of which is of multiplicity 2. The condition (2.1) excludes the root Ci = —1 
from (3.3). The remain is one root of multiplicity 2: 

(3.4) 

The asymptotic expansion (2.3) corresponding to (3.4) is 


U{p,q) =P^ 


1 + ^ Pis 


(3.5) 


The case = 1. This case corresponds to the short slant segment in the upper 
boundary of the Newton polygon in Fig. 2.2. It comprises two nodes A 410 and 
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^6 8- Therefore, substituting the expansion (2.3) with = 1 into the equation 
(1.15), we get the following equation for Cc 

^6 8 Q® + A4 10 Q^ = 0. (3.6) 

The common divisor Ci^ can be factored out from the equation (3.6). Since Ci ^ 0, 
we can remove this common divisor. Then the equation takes the form 


(3.7) 

(3.8) 


^6 8 Ci^ + ^4 10 — 0. 

Taking into account (2.18), the equation (3.7) is transformed to 

^ 0 . 

The quadratic equation (3.8) has two simple root 

Ci =p, Ci = -p, (3.9) 

The condition (2.1) excludes the root Ci = —p from (3.9). Therefore as a remain 
we have only one root, which is of multiplicity 1: 

C^=p. (3.10) 

The asymptotic expansion (2.3) corresponding to (3.10) is 


ti{p,q) =pq(i + '^l3isq '’V 


(3.11) 


The case ai = 2. This case corresponds to the long slant segment in the upper 
boundary of the Newton polygon in Fig. 3.2. It comprises three nodes Teg, ^8 4 , 
and Tio Q. Therefore, substituting the expansion (2.3) with ai = 2 into the equation 
(1.15), we get the following equation for Cp. 

Tioo a^° + A8 4 a® + T6 8 C,® = 0. (3.12) 

The common divisor Ci® is factored out from the equation (3.12). Since Ci ^ 0, we 
can remove this common divisor. Then the equation takes the form 

^10 0 Ci^ + T84 Ci^+ Aq 8 = 0. (3.13) 

Taking into account (2.18), the equation (3.13) is transformed to 

aV6C.2 + l = 0. (3.14) 

The quartic equation (3.14) has four roots. All of them are complex: 

C, = (^2 + 1)1, C^ = {V2-l)i, (3.15) 

C, = -(\/2 + l)i, C^ = -(V2-l)i. (3.16) 

Here i = \/—1. The roots (3.16) are excluded by the condition (2.1). The remain is 
two root (3.15) of multiplicity 1. They yield the following asymptotic expansions: 


U{p,q) = {V2 + l)iq'^(l + '^Pisq 

^ S^l ^ 

/ oo 

U{p,q) = {V2- l)iq'^(l + '^P^sq' 




(3.17) 


The results (3.5), (3.11), (3.17) are summed up in the following theorem. 
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Theorem 3.1. For sufficiently large positive values of the parameter q, i. e. for 
q > qmin, the tenth-degree equation (1-15) has five roots of multiplicity 1 satisfying 
the condition (2.1). Three of them ti, t 2 , and t^ are real roots. Their asymptotics 
as q ^ +00 are given by the formulas 

h'-^p'^, hr^pq. (3.18) 

The rest two roots t^ and t^ of the equation (1-15) are complex. Their asymptotics 
as q ^ +00 are given by the formulas 

t4-^ {V2 + l)[q^, t5~(^-l)ig2. (3.19) 


The complex roots (3.19) do not provide perfect cuboids. However, below they 
are important for determining the exact number of real roots. 

4 . Asymptotic estimates for real roots. 

According to the formula (3.18) the roots ti and t 2 are not growing as q ^ +oo. 
For this reason we do not need to calculate (dis in (3.5) for them. But we need to 
find estimates for remainder terms i?i and R 2 in the formulas 


ti + i?i(p,(?), t2=p'^ + R2{p,q) (4.1) 


as q ^ + 00 . Our goal is to obtain estimates of the form 


\RiiPjq)\ < , where i = 1 , 2 . 

q 

In order to get such estimates we substitute 



(4.2) 


(4.3) 


into the equation (1.15). Then we perform another substitution into the equation 
obtained as a result of substituting (4.3) into (1.15): 


(4.4) 


Upon two substitutions (4.3) and (4.4) and upon removing denominators the equa¬ 
tion (1.15) is written as a polynomial equation in the new variables c and z. It is 
a peculiarity of this equation that it can be written as 


-I- f{c,p, z) = 4 p® c^. 


(4.5) 


Here f{c,p, z) is a polynomial given by an explicit formula. The formula for /(c,p, z) 
is rather huge. Therefore it is placed to the ancillary file strategy_formulas.txt 
in a machine-readable form. 

Let (7 ^ 59p and let the parameter c run over the interval from —5p^ to 0: 

-5p^<c<0. ( 4 . 6 ) 
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From q ^ 59p and from (4.4) we derive the estimate \z\ ^ l/59p“^. Using this 
estimate and using the inequalities (4.6), by means of direct calculations one can 
derive the following estimate for the modulus of the function /(c,p, z): 

\f{c,p,z)\ <15p^^. (4.7) 

For fixed p and z the estimate (4.7) means that the left hand side of the equation 

(4.5) is a continuous function of c taking the values within the range from 

to while c is in the interval (4.6). The right hand side of (4.5) is also 

a continuous function of c. It decreases from 100 to 0 in the interval (4.6). 
Therefore somewhere in the interval (4.6) there is at least one root of the equation 

(4.5) . 

The parameter c is related to the initial variable t by means of the formula (4.3). 
The inequalities (4.5) for c imply the following inequalities for t: 

p2 _ ^P_ 

q 

The inequalities (4.8) and the above considerations prove the following theorem. 

Theorem 4.1. For each q ^ 59 p there is at least one real root of the equation 
(1.15) satisfying the inequalities (4.8). 

The above considerations can be repeated for the case where the parameter c 
runs over the interval from 0 to 5p^. In this case due to (4.3) from 

0 < c < 5p^ 


we derive the inequalities 

p^ < t < p^ + (4.9) 

9 

for the variable t and hence we obtain the following theorem. 

Theorem 4.2. For each q ^ 59 p there is at least one real root of the equation 
(1.15) satisfying the inequalities (4.9). 

Now let’s proceed to the growing root ts of the equation (1.15) (see Theorem 3.1). 
Upon refining the asymptotic formula (3.18) for t^ looks like 

16 

t3=pq - \-R3{p,q)- (4.10) 

q 

The formula (4.10) is in agreement with the expansion (3.11). It means that 

P 31 = 0, /332 = —16p^. 

Like in (4.2), our goal here is to obtain estimates of the form 

\R3{q)\ < ®. 


( 4 . 11 ) 
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In order to get such estimates we substitute 


t=pq- 


16p^ c 

- 1 - 2 

q q^ 


(4.12) 


into the equation (1.15). Immediately after that we perform the substitution (4.4) 
into the equation obtained by substituting (4.12) into (1.15). As a result of two 
substitutions (4.12) and (4.4) upon eliminating denominators the equation (1.15) 
is written as a polynomial equation in the new variables c and z. It looks like 


p{c,p,z) = -2p^c. 


(4.13) 


Here ip{c,p^ z) is a polynomial of three variables. The explicit formula for (p{c,p, z) 
is rather huge. Therefore it is placed to the ancillary file strategy_formulas.txt 
in a machine-readable form. 

Let q ^ 59p and let the parameter c run over the interval from —5p^ to 5p^: 

-5p'‘<c<5p^. (4.14) 

From q ^ 59p and from (4.4) we derive the estimate \z\ ^ l/59p“^. Using this 
estimate and using the inequalities (4.14), by means of direct calculations one can 
derive the following estimate for the modulus of the function (p{c,p, z): 

\(p{c, z)\ < 7p^. (4.15) 


For fixed p and z the estimate (4.15) means that the left hand side of the equation 

(4.13) is a continuous function of c taking the values within the range from — 

to while c is in the interval (4.14). The right hand side of the equation (4.13) 
is also a continuous function of c. It decreases from lOp® to —lOp® in the interval 

(4.14) . Therefore somewhere in the interval (4.14) there is at least one root of the 
polynomial equation (4.13). 

The parameter c is related to the initial variable t by means of the formula (4.12). 
The inequalities (4.14) for c imply the following inequalities for t: 


pq- 


16 p^ 

q 


5p^ 

<t<pq 


16 p^ 

q 



(4.16) 


The inequalities (4.16) and the above considerations prove the following theorem. 

Theorem 4.3. For each q ^ 59 p there is at least one real root of the equation 
(1.15) satisfying the inequalities (4.16). 

Theorems 4.1, 4.2, and 4.3 solve the problem of obtaining estimates of the form 
(4.2) and (4.11) for the remainder terms in the refined asymptotic expansions (4.1) 
and (4.10) for q ^ 59p. 


5 . Asymptotic estimates for complex roots. 

Let’s proceed to complex roots of the equation (1.15). Upon refining the asymp¬ 
totic formula (3.19) for the complex root t^ is written as 

t 4 = (\/2 -I- 1) iq^ -I- (v^ — 2) ip^ -I- i? 4 (p, g), where i = \/^. (5.1) 

The formula (5.1) is in agreement with the first expansion (3.17). It means that 
/34 i = 0 and /342 = (4 — 3-\/2)p^- Like in the formula (4.2) and in the formula 
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(4.11), our goal here is to obtain estimates of the form 

\R4{p,q)\<^. (5.2) 

In order to get such estimates we substitute 

t = {V2 + l)iq^ + {V2-2)ip^+'^ (5.3) 

into the equation (1.15). Immediately after that we perform the substitution (4.4) 
into the equation obtained by substituting (5.3) into (1.15). As a result of two 
substitutions (5.3) and (4.4) upon eliminating denominators the equation (1.15) is 
written as a polynomial equation in the new variables c and z. It looks like 

tl}{z,p,c) = 16 c. (5.4) 

Here tjj(c,p, z) is a polynomial of three variables with purely real coefficients. The 
explicit formula for '\j){c,p, z) is rather huge. Therefore it is placed to the ancillary 
file strategy_formulas.txt in a machine-readable form. 

Let q ^ 59p and let the parameter c run over the interval from —5p^ to 5p^: 

—5p^ < c < 5p^. (5.5) 

From q ^ 59p and from (4.4) we derive the estimate \z\ ^ l/59p“^. Using this 
estimate and using the inequalities (5.5), by means of direct calculations one can 
derive the following estimate for the modulus of the function 'ijj{c,p,z): 

\i’{c,P,z)\ <15p^. (5.6) 

For fixed p and z the estimate (5.6) means that the left hand side of the equation 
(5.4) is a continuous function of c taking its values within the range from —15p^ 
to 15p^ while c runs over the interval (5.5). The right hand side of the equation 
(5.4) is also a continuous function of c. It increases from —8Qp^ to 80p^ in the 
interval (5.5). Therefore somewhere in the interval (5.5) there is at least one root 
of the polynomial equation (5.4). 

The parameter c is related to the initial variable t by means of the formula (5.3). 
Therefore the inequalities (5.5) for c imply the following inequalities for t: 

{V2 + l)q'^ + {V2-2)p‘^-^ <lmt< iV2 + l) q^+ iV2-2)p^+ ^. (5.7) 

The inequalities (5.7) and the above considerations prove the following theorem. 

Theorem 5.1. For each q ^ 59p there is at least one purely imaginary root of the 
equation (1.15) satisfying the inequalities (5.7). 

The complex root t^ is similar to the root t 4 . Upon refining the asymptotic 
formula (3.19) for the complex root t^ is written as 


ti = {V2 - l)[q'^ + {V2+ 2)ip^ + R 5 {p,q), where i = \/^. 


(5.8) 
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The formula (5.8) is in agreement with the second expansion (3.17). It means that 
/Ssi = 0 and j3^2 = (4 + 3 \/2)p'^. Like in the formulas (4.2), (4.11), and (4.2), our 
goal here is to obtain estimates of the form 

\R5{p,q)\ < (5.9) 

q 

In order to get such estimates we substitute 


t — (v^ — 1) i + (v^ + 2)ip^H—- (5.10) 

into the equation (1.15). Immediately after that we perform the substitution (4.4) 
into the equation obtained by substituting (5.10) into (1.15). As a result of two 
substitutions (5.10) and (4.4) upon eliminating denominators the equation (1.15) 
is written as a polynomial equation in the new variables c and z. It looks like 

rj{z, c) = 16 c. (5.11) 

Here r](c,p, z) is a polynomial of three variables. The explicit formula for rj{c,p, z) 
is rather huge. Therefore it is placed to the ancillary file strategy_formulas.txt 
in a machine-readable form. 

Let <7 ^ 59p and let the parameter c run over the interval from —5p^ to 5p^ 
(see (5.5)). From q ^ 59p and from (4.4) we derive the estimate \z\ ^ l/59p“^. 
Using this estimate and using the inequalities (5.5), by means of direct calculations 
one can derive the following estimate for the modulus of the function 77 ( 0 , p,z): 

l?7(c,P,2)l < 15p^ (5.12) 

For fixed p and z the estimate (5.12) means that the left hand side of the equation 

(5.11) is a continuous function of c taking its values within the range from —15p^ 
to 15while c runs over the interval (5.5). The right hand side of the equation 

(5.12) is also a continuous function of c. It increases from —80p^ to 80p^ in the 
interval (5.5). Therefore somewhere in the interval (5.5) there is at least one root 
of the polynomial equation (5.12). 

The parameter c is related to the initial variable t by means of the formula (5.10). 
Therefore the inequalities (5.5) for c imply the following inequalities for t: 

{V2-l)q'^ + {y/2 + 2)p^-^ <lmt< {V2-1) q^+ {V2 + 2) p'^+ ^. (5.13) 

q q 

The inequalities (5.13) and the above considerations prove the following theorem. 

Theorem 5.2. For each q ^ 59p there is at least one purely imaginary root of the 
equation (1.15) satisfying the inequalities (5.13). 

Theorems 5.1 and 5.2 solve the problem of obtaining estimates of the form (5.2) 
and (5.9) for the remainder terms in the refined asymptotic expansions (5.1) and 
(5.8) for q ^ 59p. Along with Theorems 4.1, 4.2, and 4.3, they separate the roots 
^ 2 , ^ 3 , ti, ^5 of the equation (1.15) from each other for sufficiently large q and 
provide rather precise intervals for their location. 
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6. Non-intersection of asymptotic intervals. 

Theorems 4.1, 4.2, 4.3, 5.1, 5.2 define five asymptotic intervals (4.8), (4.9), (4.16), 
(5.7), and (5.13) for q ^ 59p. It is easy to see that the intervals (4.8) and (4.9) 
do not intersect. For the other pairs of intervals among (4.8), (4.9), (4.16), (5.7), 
(5.13) this is not so obvious. Therefore we need some elementary lemmas. 

Lemma 6.1. For q ^ 59p the asymptotic intervals (4.8), (4.9), (4.16), (5.7), and 
(5.13) do not comprise the origin. 

Proof. Indeed, from q ^ 59p for the left endpoint of the interval (4.8) we derive 


P^- 


5p^ 


^ 2 


54p^ 

59 


> 0 . 


( 6 . 1 ) 


The left endpoint of the interval (4.9) is obviously positive: p^ > 0. In the case of 
the interval (4.16) from q ^ 59 p we derive 


pq- 


16p3 5p4 204430p 




3481 


> 58p^ > 0. 


( 6 . 2 ) 


In the case of the imaginary intervals (5.7) and (5.13) from q ^ 59p we derive 


5 

{V2 + 1) q^ + - 2)p'^ - ^ > 8403p2 > 0, 

5 

{V2-l)q^ + {V2 + 2)p^ ~q^ 1445 ^^ > 0 . 


(6.3) 


The above inequalities (6.1), (6.2), and (6.3) prove Lemma 6.1. □ 

Lemma 6.1 means that for q 59 p the real intervals (4.8), (4.9), and (4.16) 
do not intersect with the purely imaginary intervals (5.7) and (5.13). Moreover, 
the inequalities (6.1), (6.2), and (6.3) show that all of these intervals are located 
within positive half-lines of the real and imaginary axes. Therefore any roots of the 
equation (1.15) enclosed within these intervals satisfy the condition (2.1). 

Lemma 6.2. For q ^ 59p the real asymptotic intervals (4.8), (4.9), and (4.16) do 
not intersect with each other. 

Proof. The open intervals (4.8) and (4.9) are adjacent. They have one common 
endpoint t = p"^, but this endpoint does not belong to them. Therefore the intervals 
(4.8) and (4.9) do not intersect with each other. 

In order to prove Lemma 6.2 it is sulficient to compare the right endpoint of the 
interval (4.9) with the left endpoint of the interval (4.16). From q ^ 59p we derive 


5p^ 64p^ 2 

q 59 


(6.4) 


Comparing (6.4) with (6.2), we conclude that 


5p3 IQpS 

— <pq - 

q q 


5p^ 


q 


(6.5) 
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for q ^ 59p. The inequality (6.5) completes the proof of Lemma 6.2. □ 

Lemma 6.3. For q ^ 59p the imaginary asymptotic intervals (5.7) and (5.13) do 
not intersect with each other. 

Proof. In order to prove Lemma 6.3 it is sufficient to compare the bottom endpoint 
of the interval (5.7) with the top endpoint of the interval (5.13): 

Ini(e:L - CrO =2q^-4p^- (6.6) 

From q ^ 59 p and from (6.6) we derive the inequalities 

Ini(e:L - C) > > 6957p2 > 0. (6.7) 

The inequalities (6.7) complete the proof of Lemma 6.3. □ 

Lemmas 6.1, 6.2, and 6.3 are summed up in the following theorem. 

Theorem 6.1. For q ^ 59p five roots ti, t 2 , t^, t^, t^ of the equation (1.15) obeying 
the condition (2.1) are simple. They are located within five disjoint intervals (4.8), 

(4.9) , (4.16), (5.7), (5.13), one per each interval. 

Due to (2.2) Theorem 6.1 locates all of the ten roots of the equation (1.15). 

7. Integer points of asymptotic intervals. 

It is easy to see that the asymptotic intervals (4.8), (4.9), (4.16) become more 
and more narrow if p is fixed and q —>■ +oo. Using this observation, one can easily 
prove the following two theorems. 

Theorem 7.1. If q 59p and q > 5p^, then the asymptotic intervals (4.8) and 

(4.9) have no integer points. 

Theorem 7.2. If q^ 59p and q^ > lOp^, then the asymptotic interval (4.16) has 
at most one integer point. 

The next theorem is more complicated. 

Theorem 7.3. If q 59p and q ^ 16p^ + 5p/16, then the asymptotic interval 
(4.16) has no integer points. 

Proof. Note that (1 + x)^ > 1 + 2x for any positive x. This inequality yields 

1 + a: > Vl + 2a; for any a; > 0. (7.1) 

Let’s write the inequality q ^ 16p^ + 5p/16 in the following way: 

,-8p=>8pM^ = 8p=(l + j5^). (7.2) 

Setting X = 5/(128p^) in (7.1) and applying it to (7.2), we get 

5 


q - 8p^ > 8p^ 


64p2 


(7.3) 
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Both sides of the inequality (7.3) are positive. Squaring them, we obtain 


(7.4) 


(,-8pV>(i4p»(l + ^). 

Expanding both sides of the inequality (7.4), we bring it to 

— 16p^ q > 5p^. (7.5) 

And finally, dividing both sides of the inequality (7.5) by q^, we write it as 

16p^ 5p^ 


< 1 . 


q q- 

Apart from q ^ 16p^ + 5p/16, we have the inequality q ^ 59p that yields 

A < 

q? 


(7.6) 


59 pg 59 g 


5p3 ^ I6p3 


(7.7) 


Applying (7.6) and (7.7) to (4.16), we derive the following inequalities: 

pq — 1 < t < pq. 


(7.8) 


Since pg is integer, the inequalities (7.8) have no integer solutions for t. This means 
that the interval (4.16) has no integer points. Theorem 7.3 is proved. □ 


8. A STRATEGY FOR NUMERIC SEARCH. 

The numeric search for perfect cuboids in the case of the second cuboid conjecture 
(see Conjecture 1.1) is based on the equation (1.15). The equation (1.15) is related 
to the equation (1.8) through the substitutions (1.10) and (1.13). Substituting 
either (1.10) or (1.13) into the inequalities t > a,t > b, and t > u from Theorem 1.2, 
we get the same result expressed by the inequalities 


t>p'^, t>pq, t>q'^. (8.1) 

Similarly, substituting either (1.10) or (1.13) into the inequality (a + t) (b + t) > 2t^, 
we get the same result expressed by the inequality 

(p^ + f) (pg + t) > 2t^. (8.2) 

Theorem 1.2 specified for the case of second cuboid conjecture (see Conjecture 1.1) 
is formulated in the following way. 

Theorem 8.1. A triple of integer numbers p, q, and t satisfying the equation 
(1.15) and such that p ^ q are coprime provides a perfect cuboid of and only if the 
inequalities (8.1) and (8.2) are fulfilled. 

Generally speaking, the numeric search based on Theorem 8.1 is a three-para- 
metric search. The inequalities (8.1) set lower bounds for t, but they do not restrict 
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i to a finite set of values. The inequality (8.2) is different. Since p and q are positive, 
the inequality (8.2) can be written as follows: 


t < 


p^ +pq _^p yjp'^ + &pq + 


(8.3) 


Due to (8.1) and (8.3) for each fixed p and fixed q one should iterate only over a 
finite set of integer values of f, i. e. the search is two-parametric in effect. 

Theorems 7.1 and 7.3 strengthen the restrictions. They say that for each fixed 
p one should iterate over a finite set of values of q and t, i. e. the search for perfect 
cuboids becomes effectively one-parametric. 

Theorems 7.1 and 7.3 can be further strengthened with the use of the inequalities 
(8.1). Assume that the condition q ^ 59p is fulfilled and assume that t belongs 
to the first asymptotic interval (4.8) (see Theorem 6.1). Then from (8.1) and (4.8) 
we derive two contradictory inequalities t > p^ and t < p^. Hence we have the 
following theorem. 

Theorem 8.2. If q ^ 59p, then the asymptotic interval (4.8) has no points satis¬ 
fying the inequalities (8.1). 

Similarly, if g ^ 59p and if t belongs to the second asymptotic interval (4.9), 
then from the inequalities (4.9) and (8.1) we derive 


t>q^^ (59p)2 = 3481p^ 
t<p^ + —. 

q 


The inequalities (8.4) mean that 


Q/101 2 , 2 I 

3481 p <p -I- 


(8.4) 


(8.5) 


From (8.5) one easily derives the inequality q < 5p/3480 = p/696 that contradicts 
q 1)9 p. Hence we have the following theorem. 

Theorem 8.3. If q ^ 59p, then the asymptotic interval (4.9) has no points satis¬ 
fying the inequalities (8.1). 

Finally, assume that g ^ 59p and let t belong to the third asymptotic interval 
(4.16). We know that g ^ 59p implies (7.7). From (7.7) and (4.16) we derive 


16p3 5p4 

t<pq -^- j- <pq. 

q q^ 


( 8 . 6 ) 


The inequalities (8.6) contradict the inequality t > pq from (8.1). This contradic¬ 
tion proves the following theorem. 

Theorem 8.4. If q ^ 59p, then the asymptotic interval (4.16) has no points sat¬ 
isfying the inequalities (8.1). 

The fourth and fifth asymptotic intervals (5.7) and (5.13) are purely imaginary. 
They do not provide perfect cuboids. Therefore Theorems 8.1, 8.2, 8.3, and 8.4 are 
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summarized in the following theorem. 

Theorem 8.5. If q ^ 59p, then the Diophantine equation (1-15) has no solutions 
providing perfect cuboids. 

Theorem 8.5 is a background for an optimized strategy of numeric search for 
perfect cuboids. It says that only for q < 59 p perfect cuboids are expected. 
Applying the inequality q < 59 p to the formula (8.3) we can simplify it as 

t<61p2. (8.7) 

The upper bound (8.7) is more simple than (8.3), though it can be computationally 
less time-efficient for large values of p. Along with (8.1) and the inequality 

q < 59p, (8.8) 

it provides the following very simple computer code for our optimized strategy: 

for p from 1 by 1 to -|-oo do 
for q from 1 by 1 to 59*p-l do 
if pOq and gcd(p,q) = l 
then 

for t from max(p~2,p*q,q~2)+l by 1 to 61*p~2-l do 
if Q_pqt=0 and (p~2+t)*(p*q+t)>2*t~2 
then 

Str: =sprintf ("Cuboid is found: p=7oa, q=yoa, t=yoa. " ,p,q,t) : 
writeline(default,Str): 
end if: 
end do: 
end if: 
end do: 
end do: 

Here gcd(p,q) stands for the greatest common divisor of p and q, while Q_pqt is a 
computer version of the formula (1.12). In practice the infinity sign -|-oo should be 
replaced by some particular positive integer. 

9. Concluding remarks and acknowledgments. 

Though they look very simple, Theorem 8.5 and the inequality (8.8) constitute 
the main result of the present work. As for the above code, it should be further 
optimized e. g. by some tricky algorithms for fast computing the values of Q_pqt. 
One of such further optimized versions of this code has been run on a desktop PC 
for p from 1 to 100. It took about 7 hours for that. No perfect cuboids were found. 

The authors are grateful to Mr. Seth Kitchen from Missouri University of Science 
and Technology whose e-mail letter asking about a cuboid coding idea became an 
impetus for this work. 
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